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Abstract: Several Cardassian universe models including the original, modified polytropic and 
exponential Cardassian models are constrained by the latest Constitution Type la supernova data, 
the position of the first acoustic peak of CMB from the five years WMAP data and the size of 
baryonic acoustic oscillation peak from the SDSS data. Both the spatial flat and curved universe 
are studied, and we also take account of the possible bulk viscosity of the matter fluid in the flat 
universe case. 
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1. Introduction 

Nowadays, there are many dark energy models and modified gravity theories proposed to explain 
the current accelerating expansion of the universe, which has been confirmed by the observations 
like Type la supernovae (SNe la), CMB and SDSS et al. The dark energy models assume the 
existence of an energy component with negative pressure in the universe, and it dominates and 
accelerates the universe at late times. The cosmological constant seems the best candidate of dark 
energy, but it suffers the fine tuning problem and coincidence problem, and it may even have the 
age problem [Q. To alleviate these problems, many dynamic dark energy models were proposed. 
However, people still do not know what is dark energy. 

Since the Einstein general gravity theory has not been checked in a very large scale, then one 
does not know whether this gravity theory is suitable or not for studying the observational data 
like SNe la, and maybe the accelerating expansion of universe is due to the gravity theory that 
differs from the general gravity. Thus, many modified gravity theories like f(R), DGP et al. are 
proposed to explain the accelerating phenomenology. The Cardassian model is a kind of model in 
which the Fridemann equation is modified by the introduction of an additional nonlinear term of 
energy density, and we will briefly review on this model in the next section. 

Dissipative processes in the universe including bulk viscosity, shear viscosity and heat transport 
have been conscientiously studied The general theory of dissipation in relativistic imperfect 
fluid was put on a firm foundation by Eckart||, and, in a somewhat different formulation, by 
Landau and LifshitzfQ. This is only the first order deviation from equilibrium and may has a 
causality problem, the full causal theory was developed by Isreal and Stewart H, and has also been 
studied in the evolution of the early universe]?]]. However, the character of the evolution equation 
is very complicated in the full causal theory. Fortunately, once the phenomena are quasi-stationary, 
namely slowly varying on space and time scale characterized by the mean free path and the mean 
collision time of the fluid particles, the conventional theory is still valid. In the case of isotropic 
and homogeneous universe, the dissipative process can be modeled as a bulk viscosity £ within a 
thcrmodynamical approach, while the shear viscosity rj can be neglected, which is consistent with 
the usual practice||. For works on viscous dark energy models, see ref.Q. 

The bulk viscosity introduces dissipation by only redefining the effective pressure, p e ff, accord- 
ing to p e f f — p~3£H where £ is the bulk viscosity coefficient and H is the Hubble parameter. The 
condition £ > guaranties a positive entropy production, consequently, no violation of the second 
law of the thermodynamics ffO|| . The case £ = tH, implying the bulk viscosity is proportional to 
the fluid's velocity vector, is physically natural, and has been considered earlier in a astrophysical 
context, see the review article of Gr0n|p"lf. 
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In this paper, we will focus on several Cardassian models including the original, modified poly- 
tropic and exponential Cardassian models and constrain their parameters by the latest Constitution 
Type la supernova data (SNela), the position of the first acoustic peak of the cosmic microwave 
background (CMB) from the five years WMAP data and the size of baryonic acoustic oscillation 
(BAO) peak from the SDSS data. We have consider the case of spatial flat and curved universe and 
the case of flat universe with the bulk viscosity. After a lengthy numerical calculation, we obtain 
the best fit values of the parameters in each Cardassian model. 

This paper is organized as follows: In Section 2, we present a brief review of Cardassian models, 
and derive the Hubble parameter in terms of the redshift and some parameters for several models. 
In Section 3, we analysis each model with statistical method and constrain their parameters with 
the observational data. In the last section, we give some conclusions and discussions. 

2. The Cardassian model 

Assuming the universe is homogeneous and isotropic, i.e. 

ds 2 = -dt 2 + a(t) 2 (y^2 + r 2 (d6 2 + sin 2 6dcf, 2 )^ , (2.1) 

where k is the spatial curvature, the modified Friedmann equation for Cardassian model is given 
by 

9{P) 
a 2 3 



H 2 + 4 = m , (2.2) 



where p is the total energy density of matter and radiation and we will neglect the contribution 
of radiation for the late-time evolution of the universe. In eq.( |2.2| ), the function g(p m ) reduces to 
p m in the early universe, then eq. fl2.2|) reduces to the ordinary Friedmann equation during early 
epochs such as primordial nucleosynthesis. However, it differs from the FRW universe at the redshift 
z < 0(1), during which it will gives rise to accelerated expansion. Different forms of the function 
g{Pm) corresponds to different Cardassian models, and we will focus on the original Cardassian 
model (OC) [jl2| , the modified polytropic Cardassian model (MPC) the exponential model 
(EC) Q, their flat versions ( FOC, FMPC, FEC ), in which the spatial curvature is neglected, 
and their viscous versions ( VOC, VMPC, VEC ) p^| , in which the bulk viscosity of the matter is 
taken account and the spatial curvature is also neglected. We summaries these models in Table [|. 
Recent works on constraining the Cardassian universe, see ref. uM . 
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Table 1: Summary of Cardassian models with different functions of g(p m ). Here, p car d is a character 
energy density. 
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Energy conservation of pressureless matter is given by 



3H( Pm - = , 



(2.3) 



where £ TO is the bulk viscosity for the matter p m . Following 17 , the function g could be written 
as g = p rn + pk, where pk is so called Cardassian term, which may indicate that our observable 
universe as 3 + 1 dimensional brane in extra dimensions. Thus, the total energy density can be 
written as 

g + 3H{g+pT-HH) = 0, (2.4) 

where £ is the bulk viscosity for the total energy density g(p m )- Here, px is defined as the effective 
pressure of total fluid without bulk viscosity, and the first law of thermodynamics in an adiabatic 
expanding universe gives 

dg 



PT = P-n 



dp r , 



(2.5) 



Therefore, one can get £ = -§^-£,m from eqs. (2.3) and (2.4). In the following, we will choose 
£ = tH, in which the cosmological dynamics can be analytically solvable Jl5| and r is a constant. 
Then, the conservation law (2^4) becomes 



0, 



(2.6) 



where / = g/p m , V — hi(Pm / Pcard) the prime denotes the derivative with respect to x = ln(a) 
— ln(l + z), and z is the redshift. 

For the VOC model, / = 1 + e^™ -1 ^, and the solution is 



Prn =p m0 (l + z) 3(1 ~ T) 



1 - T + (n - T) (prno/ Pcard) 



n-l-i 



1 - T + (n - T) {p m / Pcard) 



n-1 



(2.7) 



where p m Q is the present value of the matter's energy density, and the Hubble parameter E 2 
H 2 /Hq, is given by 



E 2 = fi m0 (l + z) 3 ^ 



n + i 



h + F voc (z) 



[1 + F voc (z)(^ -1)] 



where tt m o = p m o / (SH 2 ) , Hq is the present value of the Hubble parameter and 



n 



1 — QmO 



1 - T 



Here the function F voc {z) — (p m / p m o) n 1 satisfies 

(1 + 2 )3(l-r)(n-l) 



TP 



Tl + 1 

fi + F vo 



(2.8) 



(2.9) 



(2-10) 



from which one can get the solution for F voc and substitute it into eq.(2.8), then one obtains the 
Hubble parameter in terms of z and parameters Q m o,n, r. When r = 0, the solution is rather 
simple, and the Hubble parameter ( |2.8|) becomes 

e 2 = a m0 (i + zf + (i - n m0 )(i + zf n . (2.ii) 

For the VMPC model, / = (1 + e^™- 1 ^) 1 /?, and the solution is 

'l-T+(n-T)(p m o/pcard) qin - 1} 



Pm =p m0 (l + z) 3(1 T) 



l-T+(n-T)(p m /p card )i( n - 1 ) 



(2.12) 
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and the Hubble parameter is given by 



E 2 = Q m0 (l + zf^ 



r 2 + 1 



7~2 H~ F vrri p C [z^) 



[l + F„ mpc (#J-f , (2.13) 



where 



T2 = 



lL mO 
l-tt q 



1 - T 



Here the function F vmpc (z) = (p m / p m o) q ^ n ^ satisfies 



f 2 + 1 



T2+F, 



vmpc 



F — ft 4- 7 \3q(l-T){n-l) 

* vmpc ~ V ' / 

and when r = 0, the Hubble parameter ( [2.18 ) becomes 

E 2 = Q m0 (l + zf [l + (Q m q - 1)(1 + z) 3 ^"" 1 ) 
For the VEC model, / = exp (e~ ny ), and the solution is 

n - (p m o/Pcard) n (l - T) 



Prn = Pmo(l + zf 

and the Hubble parameter is given by 
E 2 = fi m0 (l + zf 

where 



n ~ (Prn/Pcard) n (l ~ t) 
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Here the function J 7 !„ ec (z) = (p m /pmo) n satisfies 



F„ ec = (l + z) 
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(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



and when r = 0, the Hubble parameter ( 2.18 ) becomes 

E 2 = fl m0 (l + zf exp ( - (1 + z)- 3 " lnO m0 
We summarize all the solutions of Hubble parameters for each model in Table |l|. 

3. Statistical analysis with the observational data 

In general, the expansion history of the universe H(z) or E(z) can be given by a specific cosmological 
model or by assuming an arbitrary ansatz, which may be not physically motivated but just designed 
to give a good fit to the data for the luminosity distance dz or the 'Hubble-constant free' luminosity 
distance Dl defined by 

H d L 



Dr = 



(3.1) 



where the light speed c is recovered to show that Dl is dimensionless. In the following, we will 
take the first strategy that assuming the Hubble parameter H(z; a±, ■ • • , a n ) with some parameters 
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(oi, • • • ,a n ) predicted by the class of Cardassian models could be used to describe the universe, 
and then we obtain the predicted value of by 



VP 



k0\ 



k-0\ 



dz' 



E(z';ai, ■ ■ ■ ,a n ) 



(3.2) 



where Qko = — ^/( a o^o) anc ^ Sinn(x) = sin(x), x, sinh(a;) for respectively a spatially closed (fife < 
0), flat (flk = 0) and open (f^ > 0) universe. 

On the other hand, the apparent magnitude of the supernova is related to the corresponding 
luminosity distance by 

<Il{z) 



fi(z) = m(z) — M = 5 log 



10 



Mpc 



25, 



(3.3) 



where fi(z) is the distance modulus and M is the absolute magnitude which is assumed to be 
constant for standard candles like Type la supernovae. One can also rewrite the distance modulus 
in terms of Dl as 

fi(z) = 51og 10 D L (z) + , (3.4) 



where 



fl = 5 log! 



25 



-51og 10 



h + 42.38 , 



(3.5) 



V Mpc 

is the zero point offset, which is an additional model independent parameter. Thus, we obtain the 
predicted value of [i th by using the value of and the observational value of ix ohs we used is the 
latest data called the constitution data JL8[ , which contains 397 data points including the the 307 
Union data set and 90 CFA data set. 

There are also some constraints from CMB and BAO observations. We will take the parameter 
R from the CMB data |^TJ and the parameter A from the SDSS data [gl] as well as the supernova 
data to constrain parameters of the Cardassian models. The parameter R is defined as 



R = 



9. 
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(3.6) 



where zi s = 1090 is the redshift of the last scattering surface and the observational value is given 
by Robs — 1-170 ± 0.019. While the parameter A is defined as 



mO 
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(3.7) 



where z\ = 0.35 and the observational value is given by A obs = 0.469(0. 95/0. 98)" 35 ± 0.017. 

In order to determine the best value of parameters (with la error at least ) in the Cardassian 
models, we will use the maximum likelihood method and need to minimize the following quantity 



-2,2 i 2 

X.SN + XCMB + XBAO i 
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where <ii is the la error of the observation value n° bs (zi). Since the nuisance parameter /io is 
model-independent, then we analytically marginalize it by using a flat prior P(fM)) = 1: 



-2 In 
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(3.11) 



(3.12) 
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(3.14) 



then, from now on, we will work with Xsn an( ^ to uiinimize \ 2 m eq.(3.8). The best fit parameter 
values and the corresponding Xmin an d Xmin/ DOF will be summarized for each model in Table |^. 
Here DOF is the degree of freedom defined as 



DOF = N - 



(3.15) 



where N is the number of data points, and v is the number of free parameters. 

We now apply the maximum likelihood method for each model in Tabic |l], and we summary 
the results in Table |^ including the minimum values of \ 2 an d the best fit parameters with la 
confidence level for each model. Since both the FOC and FEC model contain two parameters, we 
also plot the contours from la to 3a confidence levels for them, see Fig. [l] and Fig. ||. For each 
model, the predicted dimensionless luminosity is plotted in Fig. |^, from which one can see that 
these Cardassian models predict almost the same luminosity distance with their best fit parameters. 
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Table 2: Result: The minimum value of \ 2 an d the best fit parameters with la confidence level in each 
model. 



4. Discussion 

We have used the 397 SNe la, CMB and SDSS data to constrain several Cardassian models. Wc 
have summarized these model in Table [l], in which different forms of the function g(p m ) have been 
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Figure 1: FOC: Constraints on Q m o and n from lcr to 3<r confidence level obtained by using 397 SNe la 
+ CMB + BAO for the FOC model. The best fit point corresponds to Q m0 = 0.270, n = 0.053. 



chosen and the corresponding Hubble parameters are also given. In particular, we discuss the 
viscous Cardassian models in Section 2., in which we rewrite the Hubble parameter in a continent 
way to do the statistical analysis. 

The fitting results are presented in Table 0, in which we have shown the minimum value of \ 2 
and the minimum value Xmin P er degree of freedom. The best fit parameters with la confidence level 
for each model are also presented in Table [2[ from which one can see that, the latest observational 
data can not distinguish these models at this classical level. In other words, they predict almost 
the same evolution history of the universe and we need to take the perturbation of universe into 
account that will be studied in our further work. 



In fact, the minimal of x 2 in eq^S^) is very sensitive to the observational error of the distance 
modulus. Once the error is smaller in the future data than that at present, not every model will fit 
the data well, then one can distinguish these models and even rule out some of them. Thus, more 
precise data are very needed. 

Since in Cardassian universe, one can expalain the accelerating expansion without introducing 
any dark energy component, it is very interesting and worth further studying. We also hope that 
future observation data could give more stringent constraints on the parameters in the Cardassian 
model. 
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